Sensitivity plots for WIMP direct detection using the annual 

modulation signature 



S. Cebrian, E. Garcia, D. Gonzalez, I. G. Irastorza, A. Morales, J. Morales, A. Ortiz de 
Solorzano, A. Peruzzi, J. Puimedon, M. L. Sarsa, S. Scopel, J. A. Villar 
Laboratorio de Fisica Nuclear. Universidad de Zaragoza 
50009, Zaragoza, SPAIN 

oY 

ON 

On . Abstract 



o 

CD 

Q 



Annual modulation due to the Earth's motion around the Sun is a well known 

signature of the expected WIMP signal induced in a solid state underground 

detector. In the present letter we discuss the prospects of this technique on 

statistical grounds, introducing annual-modulation sensitivity plots for the 

WIMP-nucleon scalar cross section for different materials and experimental 

^j- \ conditions. The highest sensitivity to modulation is found in the WIMP 

ON 

m 

fSj ' from the choice of the astrophysical parameters, while the lowest values of 

the explorable WIMP-nucleon elastic cross-sections fall in most cases within 
. one order of magnitude of the sensitivities of present direct detection WIMP 

searches. 
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U ■ I. INTRODUCTION 

Non barionic Cold Dark Matter is a basic component in cosmological 
models of structure formation. Taking into account the result of microlens- 
ing surveys, it can be stated that dark baryons can account for at most only 
about one-third of the estimated density of our dark halo. The best candi- 
dates to provide the rest are Weakly Interacting Massive Particles (WIMP), 
among which the neutralino, provided by the Supersymmetric extension of 
the Standard Model, is one of the favourites. 

It is well known that all WIMP direct searches, which look for the WIMP 
elastic scattering off the nuclei of a suitable detector, are essentially con- 
strained by the fact that the predicted differential rate of the signal has a 
dependence with decreasing energy which is hardly distinguishable from the 
background recorded in the detector. While this fact does not prevent us 
in extracting upper bounds of the WIMP-nucleus interaction cross-section 
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for each WIMP mass, a distinctive signature is needed to claim a positive 
identification of the WIMP. In any case, the non appearance of the genuine 
signature looked for could provide a background-rejection method useful to 
improve the experimental sensitivity. 

The only identification signatures of the WIMP explored up to now are 
provided by the features of the Earth's motion with respect to the Dark 
Matter halo. In particular the annual modulation effect M is provided by the 
combination of the motion of the solar system in the Galactic rest frame and 
the rotation of the earth around the Sun. Due to this effect the incoming 
WIMP's velocities in the detector rest frame change continuously during the 
year, having a maximum in summer and a minimum in winter. 

Several experiments have already searched for this effect one of which 
H reports a possible positive signal. The present situation is no doubt excit- 
ing, since experimental sensitivities of underground detectors are entering for 
the first time the supersymmetric parameter space, and a host of new exper- 
iments will soon start to probe it with even higher sensitivity. In the present 
paper we discuss on purely statistical grounds the perspectives of modulation 
searches taking into account the features of present and future experiments. 
Sensitivity plots will be introduced in order to convey in a compact way all 
the necessary information. 

II. EXTRACTING THE MODULATION SIGNAL 

The procedure to extract a modulated signal with a given period and 
phase from a set of measured count rates has been discussed by several authors 
§i§. 

Due to the Earth's rotation around the Sun, the expected count rate of 
WIMP's scatterings off the target's nucleus changes periodically in time. The 
dependence may be approximated by a cosine function with period T=l year 
and phase to = 2 nd june: 

S = S + S m cosoj(ti - t ) (1) 

where So and S m are the constant and the modulated amplitude of the signal 
respectively. The oscillating frequency is u = 2ir/T and the i index indicates 
the day. The theoretical inputs introduced in the evaluation of the WIMP- 
nucleus elastic scattering (cross sections, nuclear form factors, scalar, vectorial 
or spin-dependent nature of the coupling) and all the parameters entering in 
the halo velocity distribution (r.m.s. and escape WIMP velocity, local dark 
matter density) make the evaluation of the functions So and S m rather model 
dependent §. As is costumary, the amplitudes So and S m are expressed 
in terms of the WIMP mass mw and of the point-like WIMP-nucleus cross 
section a. 

Given a set of experimental count rates Nik representing the number of 
events collected in the i-th day and k-th energy bin (the formalism can be 



2 



easily generalized to the case of a multiple-crystal set-up), the mean value of 
N ik is: 



< N ik >= IMk = [bk + S ,k + S m> k cosuj(ti - t )] • W ik 

where the bk and the Sik = #o,ft + S m> k cosufa — to) represent the average 
background and the signal respectively, in number of counts per unit of detec- 
tor mass, time and interval of collected energy E (which is related to the recoil 
energy by the relation E = QEr where Q is the quenching factor of the 
detector). Wik = MATiAEktk are the corresponding exposures, where M is 
the mass of the detector, AEk is the amplitude of the k-th energy-bin, while 
ATi represents the i-th time bin (in the following we will assume all ATj= 
1 day). The are efficiencies that have to be taken into account whenever 
some subtraction method is used with the data. In the following they will be 
neglected. For simplicity also to w dl be omitted in the following equations. 

A general procedure to compare theory with experiment is to use the 
maximum-likelihood method. The combined-probability function of all the 
collected Nik, assuming that they have a poissonian distribution with mean 
values fiik, is given by: 

Nit 

Assuming that the average background rates bk and the efficiencies are 
constant in time, a possible WIMP oscillating solution can be searched for in 
the data. This implies that the time-dependent fluctuations of background 
and efficiencies should stand well below the size of the searched modulation 
effect, if one wishes to disregard any other less exotic explanation of a modula- 
tion of the data. If this condition is verified, a WIMP analysis is justified, and 
the most probable values of mw and a maximize L or, equivalently, minimize 
the function: 

y(mw, c) = —2 log L — const 

= 2u-2^2 Nik log [bk + S 0) k + S mt k cos uU] 

ik 

where \i = J2ik n ik and all the parts not depending on mw and a may be 
absorbed in the constant because are irrelevant for the minimization. 

The function y is minimized in a two-step procedure, first with respect to 
the time-independent parts fk = i>fc+<Sb,fc and then with respect to mw and a. 
In the first minimization the bk are free parameters with the only constraint 
to be positive. So the condition S ,fc — fjk,min> or otherwise So : k=fjk,min, is 
imposed. 



Following a standard procedure [10|, a region of n standard deviations 



around the minimum in the plane (a,mw) can be found by imposing the 
condition y(a,mw) - Vmin < n 2 . 



III. DATA REDUCTION 



Although the likelihood function L depends on all the collected data N^, 
only suitable combinations of the N^. enter in the determination of the pa- 
rameters a and mw- By expanding y in powers of x k = S m ^/ f k < few % up 
to terms of the order x\ the following expression can be obtained (here and 
in the following Var indicates the variance): 



where f3 k = ^' Wk , a k = ^' Nk , k = ^ Nk , N k = Ei N ik , 

W k = J2iWi k . The last term of Eq.([|) may be dropped, since it does not 
depend on the fitting parameters. This implies that the likelihood function 
has an approximate factorization L ~ F\{N ik ) ■ F2(X' k , a,mw) and the infor- 
mation to determine a and mw is only contained in the X' k up to negligible 
corrections. As expected on general grounds, Eq.(||) reduces asymptotically 
to a x 2 : this happens when the X' k s have gaussian distributions (in the ex- 
amples that will be discussed in the following we have checked that this is 
always verified). 

The expressions @ of the X' k s are asymptotically equal to the cosine 
projections introduced in ref. || (generalized to the case of discontinuous 
data taking and unpaired days) to extract a modulated signal, 

v _ J2i N ik cos ujU - N k (3 k 

Afc = 777-7 oT\ > < Afc > ~ ^ m ' k 



Var{X k ) = 



E N ik cos 2 ujti + !3 2 k N k - 2j3 k g N ik cos utj 
W 2 k {a k -PlY 



with a k = ^ W w° S ^ • (^ or ^ ne sa ^e of comparison with Eq.(j6|) 
it is convenient to use the identity N k = (E Ni k cos 2 Luti + $ k N k — 
2/3fc E Nik coscuti)/(a k — fl k ) in the expression of Var{X' k )). This implies 
that the likelihood function y has also the asymptotic behaviour: 

/ \ 2/ \ 2/ \ (Sm,k( CF ^ rn w) - X k ) 2 
y{cr,m w ) = x (cr,mw) ">X \ a > m W) = / . —rj — ttti ■ 

In the case of paired days {(5 k = 0) and if only time intervals close to the 
maximum and the minimum of the cosine function are considered (cos ojti ~ 
±1) the X k s reduce to the june-december differences in the collected count 
rates in each energy bin E k < E < E k + AE k . 



IV. STATISTICAL SIGNIFICANCE OF THE SIGNAL 



Once a minimum of the likelihood function is found, a positive result 
excludes absence of modulation at some confidence level probability. This 
can be checked by evaluating the quantity 5 2 = y[a = 0) — y(a, mw) m in to 
test the goodness of the null hypothesis. In order to study the distribution of 
5 2 we make use of the asymptotic behaviour (|i"o|). This implies: 



5 2 = y(p = 0) - y min ~ x 2 {<? = 0) - Xmin- 



In the case of absence of a modulation effect (i.e. time-independent Poisson- 
fluctuating data) numerical simulations show that the quantity 5 2 belongs 
asymptotically to a \ 2 distribution with two degrees of freedom. We explain 
this by the fact that once the cross section a is set to zero the likelihood 
function L no longer depends on mjy (all the So and S m functions vanish) 
and this is equivalent to fixing both the parameters of the fit at the same 
time. 

In presence of the signal the minimization procedure of y described in 
the previous section may be carried out semi-analytically using Eq.(|l0|). One 
finds for a the estimator a es t given by: 
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Var(a, 
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k Var(X k ) 
k Var{X k ) 



Var{X k ) J 

CT C st is a function of the data, such that < cr cst >= a (where brackets indicate 



mean value). In Eqs.(12,13) S mj k(mw) = S m ^(cr,mw) / & and is the 



WIMP mass that maximizes the function 5(mw), which is given by: 



v-i s m , k (mw) 2 \/Var(a est ) 

^k Var(X k ) 

From equation ( |14D 5 = <5(m^) may be interpreted, as expected, as the num- 
ber of standard deviations of the signal, with: 

Var{5) = 1; 

< 6 > =< a ^ > = V< 5 2 > -1 
y/Var(p- e st) 

so in presence of modulation 5 2 has the asymptotic distribution of a non 
central \ 2 with one degree of freedom and non centrality parameter given by 
(< a est /y/Var(a est ) >) 2 =< S 2 > -1. 
Using Eqs.(|T0|,ll) we obtain: 



S m ,k(cr, m w ) 2 A E k 



where the same days of data taking have been assumed for all the energy 
bins, and in the expression of Var{Xk) the following approximations have 
been made: 

< Nik cos 2 uti > ~ < Nik > ^ cos 2 ut i 

i i 

<N ik > ~W k (bk + S ) 

while we define the factor of merit a = A J2i cos 2 u>ti (a=l in case of a full 
period of data taking) and the terms depending on the (3k have been neglected. 

Equation (|l7|) allows to estimate the needed exposure MTa in order to 
obtain a given value of < 5 2 >. (It is worth noticing that the background bk 
indicates here the amount of counts expected from radioactive contamination 
and noise, for instance evaluated by making use of a Montecarlo simulation, 
and not the average counts, that contain also the signal). Whenever the 
distribution of 5 2 for a given WIMP model is sufficiently far apart from that 
of the background, a given experimental result may be discriminated to belong 
to one population or the other. Once a required < 5 2 > is chosen, a sensitivity 
plot may be obtained by showing the curves of constant MTa in the plane 
m-w-o. 

The statistical interpretation of the sensitivity plot is obtained from the 
degree of overlapping between the distributions of 5 2 in the two cases of ab- 
sence and presence of modulation. The situation is summarized in Table |, 
where < 5 2 > is tabulated as a function of the fraction of experiments where 
absence of modulation can be excluded (rows), and the required Confidence 
Level (columns). For instance, if < 5 2 >=14.9, there is a 90% probability to 
measure a value of 5 2 higher than 7: this range would exclude the absence of 
modulation at least at the 95% C.L. In the same way a less demanding value 
of < 5 2 >=5.6 would give a 50% probability to see an effect at the 90% C.L. 
or more. In the following section, we will adopt these two representative val- 
ues in order to discuss the prospects of modulation searches. This procedure 
allows to see in a compact way the masses and exposures needed to explore 
each different region of the WIMP myj-o parameter space. 

The expression of < 5 2 > in Eq.(|l7]) deserves a short comment. Some 
attention has been devoted in discussions on the perspectives of modulation 
searches |§Q to the fact that since the <S m ,fc parameters can vanish or be neg- 
ative, cancellations may arise in the calculation of the signal-to-noise ratio, 
so that an optimal choice should be found for the energy threshold or for the 
amplitude of the energy interval where the signal is integrated. According 



to Eq.(17), in order to exploit all the information available, the signal-to- 
noise ratio should be calculated over the whole energy spectrum, using the 
smallest bin width allowed by statistics and resolution, and then combined 
quadratically, i.e.: 
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(signal/noise)^ ombined = ^(signal/noise) 2 

i 

where i indicates the energy bin. Of course in this procedure it is crucial that 
the count rates of different bins are not correlated. 



V. SENSITIVITY PLOTS AND QUANTITATIVE DISCUSSION 

In this section we give a quantitative discussion of the sensitivity plots 
by considering different target materials: Ge, Te02, Nal with quenching 
factors, experimental thresholds and resolutions summarized in Table || (they 
have been chosen in order to be indicative of running or planned experiments 
P,|ll"|-|l4| ) . For simplicity, in order to calculate the sensitivity plots, a flat 
background will be always assumed and in the following the index k of will 
be dropped. 

As far as the astrophysical parameters entering in the calculation are 
concerned, we assume as usual that the WIMP velocities in the halo follow 
a Maxwellian distribution. The WIMP r.m.s. velocity v vms is related to the 
measured rotational velocity of the Local System at the Earth's position v\ oc 
by the relation t; 2 ms = | vf oc , while the Sun's velocity v sun in the galactic rest 
frame is given by v sun ~ (v\ oc + 12) km sec -1 (here the motion of the solar 
system with respect to the Local System is taken into account) and the Earth 
velocity is: 

dearth = V sun + V olh sin 5 COS [u(t - t Q )] 

where v or b — 30 km sec -1 and sin 5 ~ 0.51 || (5 is the angle between the 
Ecliptic and the Galactic plane). 



As it has been pointed out in the literature [15| the present uncertainty 
on v\ oc can affect the final result of WIMP direct detection calculations in 
a significant way. In the following we will vary v\ oc in its physical range, 
floe = (220 ± 50) km sec" 1 . Finally, we adopt an escape velocity v esc =650 
km/sec and take for the local halo mass density the value p=0.3 GeV/cm 3 . 

In order to compare the results for different target materials in the same 
plots, we show our results in the plane mw — a^™ 1 ' , where cxW is the WIMP 
cross section rescaled to the nucleon by adopting a scalar-type interaction 
(such as the one that would be dominant in the case of a neutralino). In the 
case of a monoatomic target this implies the transformation: 

o = a^A 2 ^ 

where A is the target atomic number, nw,N is the WIMP-nucleus reduced 
mass and nw,n the WIMP-nucleon reduced mass. Generalization to a multi- 
target species is straightforward. 

In Figures 1-5 we give an estimate of the minimal exposures needed to 
explore the WIMP parameter space by calculating the sensitivity plots for 



< 5 2 >=5.6 and i>i oc =220 km sec -1 . In all the figures the curves are obtained 
for values of MTa ranging from 10 kg- year to 100 kg- year in steps of 10 
(from top to bottom). The closed contour and the cross indicate respectively 
the 2a C.L. region singled out by the DAMA modulation search experiment 
and the minimum of the likelihood function found by the same authors ||. 
Note that by Eq.(|i~7|) each sensitivity plot corresponds to a fixed value of the 
quantity (< 5 2 > -2) /MTa, so the curves may be easily rescaled to different 
values of < 5 2 > by a suitable change in the corresponding exposures. As 
expected, the sensitivity is roughly proportional to (MT) -1 / 2 . 

Figures 1-2 refer to Ge detectors with background b=0.01 cpd/kg/keV 
and energy thresholds E t h=2 and 12 keV respectively (thresholds already ob- 
tained in the COSME [|l(| and Heidelberg-Moscow [11] experiments, respec- 
tively, and a background not far from that obtained by Heidelberg-Moscow 
and CDMS (l3|]). In Figure 3 the intermediate situation of a germanium 
detector with b=0.1 cpd/kg/keV and E^=4 keV is given (except from the 
assumption of flat background, the case depicted in Fig. 3 corresponds to 
the present performances of the IGEX experiment |l2|). The example of a 
Te02 detector is shown in Figure 4 (with b=0.01 cpd/kg/keV and E^=5 
keV, which are the foreseen performances of CUORICINO [jl4|]) while an Nal 
detector with b=0.1 cpd/kg/keV and E th =2 keV (practically the DAMA g 
values) is depicted in Figure 5. We recall once more that the values of b 
quoted for the aforementioned experiments are the levels of counts only due 
to the background, i.e. they do not include the possible signal. 

In order to be sensitive to a possible halo WIMP, sensitivity plots need to 
lie below the current most stringent upper limits to the WIMP-nucleon scalar 
cross section a^ n \ Figures 6-9 show the corresponding required minimal 
exposures as a function of the WIMP mass mw for the cases of Ge, Te02 
and Nal and calculated for the exclusion plot of Ref. [17], obtained with data 
statistically discriminated by pulse shape analysis. In each plot, dotted, solid 
and dashed curves refer to v\ oc = 170, 220 and 270 km sec -1 respectively, and 
in each case the different examples of b=0,0.01 and 0.1 cpd/kg/keV are given 
from bottom to top. 

The effect of the uncertainty on v\ oc is to displace horizontally the sensi- 
tivity plots to higher WIMP masses for lower values of v \ oc . This implies that 
for a given value of MT the upper limit on the explorable WIMP mass can 
significantly rise if v\ oc is taken in its lowest range. On the other hand, the 
curves of Figs. 6-9 do not depend on p\ oc , and are almost insensitive to Vesc: 
unless for very low WIMP masses. 

The results of the present discussion are summarized in tables D I - V\. In 
order to show the dependence of theoretical expectations on v\ oc , in the three 
tables the values t>i oc =220, 170 and 270 km sec -1 are used respectively. In 
the tables we add examples with other values of the energy thresholds and 
the background b (the complete set is given in table 0). 

In the second column of tables 1I1-V| the minimal MTa values that are 
necessary to explore the regions of the m^-a^ plane below the exclusion 
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plot of Ref. are given for < 8 >= 5.6, followed, if necessary, by a cor- 
responding interval for the WIMP mass. If no interval is shown, the given 
exposure allows to reach the exclusion plot over the whole considered range 
10 GeV< m w <1000 GeV. 

In the last two columns of tables IIIV| we summarize the prospects of 
exploration of the DAMA region. The values of MTa given in the third 
column correspond to the lowest values that give a < 5 2 >=5.6 sensitivity 
plot encompassing all the 2a contour. Finally, in the last column the same 
values are given in the case of < 5 2 >=15. 



VI. CONCLUSIONS 

In the present paper the prospects for direct searches of a WIMP mod- 
ulation effect are discussed in the case of different target materials and as a 
function of the detector mass M and the time of exposure T. 

Given its purely statistical nature, such a discussion has not taken into 
account the problems related to the systematics of a real experiment, ex- 
acerbated by the challenge of detecting a small signal depending on time. 
Consequently, although the exposure MT planned or reachable in a projected 
experiment can be an indication of its future chances, the real prospects of a 
given modulation search rely in no lesser extent in its realistic possibilities of 
keeping background stability under control for very long periods of time. For 
instance, in the examples previously discussed the signal typically amounts 
to a fraction between 1% and 5% of the average count rates, concentrated in 
the low-energy range of the spectrum. This implies that the corresponding 
WIMP models would be explorable only with a control of the stability of the 
background substantially below that range. 

We also remark that in our calculation a given sensitivity is associated to 
a value of the product MTa, irrespective to the actual number of the accu- 
mulated periods. This is due to the fact that the signal/noise ratio discussed 
in the text scales as 1/y/N with N the total amount of counts in case of 
Poissonian fluctuations (the contribution of data collected in the days near 
the maximum or the minimum of the signal is enhanced by the factor of merit 
a). However in a real experiment the collection of several periods would be 
preferred in order to have a higher control of systematics effects. 

An important feature of all the plots is that the sensitivity to modulation 
is generally a decreasing function of the WIMP mass, the highest sensitivities 
corresponding roughly to the interval 10 GeV < m^y < 130 GeV. However, 
the actual upper limit of the region of WIMP masses within the reach of a 
given experiment depends in a crucial way from the choice of the astrophysical 
parameters, something already discussed in the literature [15]. 



The optimal energy intervals where the signal should be looked for can be 
found for instance by requiring that the sensitivity plot does not change more 
than 5% over the whole range of M\y when more energy bins are added. With 
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such proviso we find 2 4- 10 keV< E < 35 keV for Ge, 2 4 5 keV< £ < 30 keV 
for Te02 and 2 keV< E < 8 keV for Nal, the lower bounds depending on the 
WIMP mass. 

The region singled out by the DAMA experiment is within the reach of 
many realistic set-ups. This can be quantified by looking at the Figs. 1-5 (for 
«ioc=220 km sec -1 ) and is shown in tables |III^[V| (also for -ui oc =170 km sec -1 
and 270 km sec -1 ). For instance, for v\ oc =220 km sec -1 , the typical needed 
exposure for a germanium detector turns up to be MTa ~ 50 kg-year in the 
two realistic examples of E t ^=2-A keV, b=0.1 cpd/kg/keV and -E t h=12 keV, 
b=0.01 cpd/kg/keV, and lowers to about MTa ~ 15 kg-year in the more 
optimistic case of i? t h=2-4 keV and b=0.01 cpd/kg/keV. 

However, as can be seen from Tables [IV] and [V], for all the considered 
nuclear targets these estimations depend in a critical way on the value of v\ oc , 
detection chances being systematically better for lower values. 

As far as TeO"2 is concerned, the achievable sensitivities seem promis- 
ing over the whole considered WIMP mass interval 10 GeV< mw is 1000 
GeV, and an exposure as low as ~ 25 kg-year could be sufficient to start 
a modulation search exploration if b=0.01 and E t /j=2 keV. However, when 
comparing these results with other kinds of detector, it is worth noticing that 
these sensitivities are very dependent on the background and threshold, and 
our assumptions are to be considered as future goals since high mass TeC>2 
bolometers are still at the R&D stage of development. 

Turning to Nal detectors, they are the first that have been used to search 
for modulation ^||, essentially due to the possibility to reach high masses. 
Their prospects for detection are a sensitive function of the obtainable thresh- 
old, since the most part of the signal is contained in the first energy bins, E < 5 
keV, because of the low quenching factor. As can be seen from Fig. 9, their 
sensitivity is higher for relatively light WIMP's, and the minimal required 
exposures turn out to be very steep functions of the WIMP mass, the heav- 
ier explorable values of mw depending critically on v\ oc and reaching up to 
m w ~ 130 GeV. 

The experimental region singled out by Ref. || (30 GeV< mw < 130 GeV 
for ^i oc =220 km sec -1 ) extends well beyond the upper limit for the detectable 
WIMP masses implied by the curves of Fig. 9 (mw ;$ 70 GeV for the same 
value of fioc)- An explanation of this fact is that the experimental region 
encompasses configurations with much lower probabilities (corresponding to a 
2-0" Confidence Level) than the ones required in the sensitivity plots. For mw 
and <j( n ) corresponding to the minimum of Ref. || and assuming MTa=50 
kg year, 6=0.5 cpd/kg/keV and E th = 2 keV, we find < 5 2 >~ 7.2. The 
experimental result published by Ref. || is 5 2 = 8.23. 

The overall picture arising from the examination of Figs. 1-5 and the 
results tabulated in Tables [Hl|-|v| turns out to be dependent, as already men- 
tioned, on the choice of the astrophysical parameter v\ oc . However a general 
trend can be drawn, in which prospects of modulation searches seem promis- 
ing provided that the WIMP signal is not far below present sensitivities. The 
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lowest values of explorable fall in most cases in the typical range of few 
xlO -10 nbarn. This value can lower up to one order of magnitude in some 
extreme cases (for instance, for MTa=1000 kg-year, a low threshold of a few 
keV and setting b=0). 

The impact that future improvements on exclusion plots would have on 
the prospects of modulation searches can be estimated in the following way. 
In Eq. ( |T7| ) , at fixed < 5 2 > , when in the denominators the inequality bk > 5b, A; 
holds, MTa is (approximately) proportional to a 2 . On the other hand, when 
bk < Sb, fc> MTa is proportional to a. This implies that an improvement of 
one order of magnitude of the exclusion plot would rise the curves of Figs. 
6-9 between one and two orders of magnitude, depending on the background 
achieved by the modulation search. 

As a last remark, we make a few comments about the case of a null 
result, when modulation can be used for background subtraction, improving so 
the exclusion plot one would obtain from time-integrated counts. Numerical 
results show that the effectiveness of this technique improves the exclusion 
plot only in case of relatively high count-rate levels. This is confirmed by 
numerical inspection of Eq. (|17D : in order for a modulation search experiment 
with background b\ to extract a better exclusion plot than another experiment 
with background b2, the following approximate condition must hold 



MT/(kg year) yEfe S m,k AE > 



< \ , f2 w^ F ( m w) ; F(mw) 



2(<5] irr1 — 2) Subthreshold 

where b\ and 62 are in cpd/kg/keV, S 2 im is the upper limit to < 5 2 > pro- 
vided by the null modulation search, the S m ,k's refer to experiment 1, and 
So, threshold refers to experiment 2, assuming that its upper limit is driven 
by its background at threshold. Typically, at the 90% C.L., one can assume 
<5ii m ~ 10, while a numerical check in experimental situations analogous to the 
ones analyzed in the previous sections gives F(mw) ^0.1 and usually much 



lower. In this case Eq.(23) implies 61 < b\ x 0.1 x MT, that, for instance, 
for b\ = bi = b and MT=100 kg year implies that a modulation analysis can 
improve the exclusion plot only if b > 0.1 cpd/kg/keV. 
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TABLES 

TABLE I. Parameter < 5 2 > as a function of the required Confidence Level (columns) and the 
fraction of successful experiments(rows). 





90% 


95% 


99% 


99.5% 


99.9% 


50% 


5.6 


7.0 


10.2 


11.6 


14.8 


60% 


6.8 


8.3 


11.8 


13.3 


16.8 


70% 


8.1 


9.9 


13.7 


15.3 


19.0 


80% 


9.9 


11.8 


16.0 


17.8 


21.8 


90% 


12.8 


14.9 


19.6 


21.6 


26.0 


95% 


15.4 


17.8 


22.9 


25.0 


29.8 


99% 


21.0 


23.8 


29.7 


32.2 


37.5 



TABLE II. Summary of the experimental parameters assumed in WIMP modulation searches 
for the different target materials discussed in the text. They are inspired to the running or projected 
experiments whose reference is shown in the last column. 





E th (keV) 


FWHM(keV) 


Quenching 


Ref. 


Ge 


2;4;12 


1 


0.25 


[11-13| 


Nal 


2 


2 


0.09 


1 


Te0 2 


2;5 


2 


0.93 


11 
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TABLE III. Summary of minimal exposures, all in kg • year, for v\ oc = 220 km sec -1 . 
indicates the energy thresholds expressed in keV, b the background (assumed constant in energy) in 
cpd/kg/keV. Exposures are estimated for the WIMP mass range in parenthesis if shown, otherwise 
over the range 10< mw ^1000. 





Exploration of not 


DAMA region: 


DAMA region: 




excluded regions (< 5 2 >= 5.6) 


< 5 2 >= 5.6 


< 5 2 >= 15 


Ge 








Ei h =2, b=0.01 


35 


15 


55 




20 (m w £110) 






E th =2, b=0.1 


80 


50 


175 




55 (m w <110) 






E th =l2, b=0.01 


25 (40< m w <110) 


50 


190 


p , 10 u n 1 

^th — J-^j D— U.l 




OOU 


1 1 on 


E th =4, b=0.1 


155 


50 


190 




65 (m w <110) 






Te0 2 








Et h =2, b=0.01 


20 


25 


90 


E th =2, b=0.1 


40 


55 


210 


Et h =5, b=0.01 


40 


40 


150 


E th =5, b=0.1 


85 


120 


435 


Nal 








E t/l =2, b=0.1 


50(m w <70) 


180 


660 
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TABLE IV. The same as in Table III for vi oc =170 km sec . 





Exploration of not 


DAM A region: 


T~\ A A /T A .... „~Z _ __ . 

JJA1V1A region: 




excluded regions (< <5 2 >= 5.6) 


< >= 5.0 


< o >= 15 


Ge 








E th =2, b=0.01 


20 


8 


30 










E th =2, b=0.1 


50 


25 


90 




20 (m^y <110) 






T — \ -f C~\ 1 A A-| 

E th =12, b=0.01 


30 (m w >45) 


25 


95 


E tfe =12, b=0.1 


1 Of / \ of \ 

135 (m w >35) 


160 


570 


&th=% b=0.1 


on 


2,0 


on 
90 




20 (mw <110j 






Te0 2 








E th =2, b=0.01 


15 


15 


50 


Et fc =2, b=0.1 


25 


30 


110 


E* h =5, b=0.01 


25 


20 


80 


E th =5, b=0.1 


50 


65 


230 


Nal 








E th =2, b=0.1 


50(m w £125) 


100 


355 



TABLE V. The same as in Table |ffl| for vi oc =270 km sec -1 





Exploration of not 


DAMA region: 


DAMA region: 




excluded regions (< S 2 >= 5.6) 


< S 2 >= 5.6 


< 5 2 >= 15 


Ge 








E th =2, b=0.01 


50 


25 


90 




40 (m w £110) 






E th =2, b=0.1 


105 


80 


285 


E th =12, b=0.01 


70 (25< m w <110) 


70 


260 


E th =12, b=0.1 


300 (25< mw <H0) 


480 


1735 


E th =A, b=0.1 


210 


85 


310 




180 (m w <H0) 






Te0 2 








E th =2, b=0.01 


30 


40 


140 


E th =2, b=0.1 


60 


90 


330 


E th =5, b=0.01 


60 


65 


230 


E th =5, b=0.1 


125 


190 


680 


Nal 








E th =2, b=0.1 


180(m w <60, m w >200) 


305 


1100 



15 



FIGURES 




M w (GeV) 

FIG. 1. Sensitivity plots in the a^-mw plane for a germanium detector with threshold energy 
E i / l =2 keV, flat background b=0.01 cpd/kg/keV and calculated for < 5 2 >=5.6. The set of curves 
correspond to different values of the exposure, MTa=10 to 100 kg-year in steps of 10 from top to 
bottom. The closed contour represents the 2a C.L. region singled out by the modulation analysis 
performed by the DAMA experiment || and the cross indicates the minimum of the likelihood 
found by the same authors. 
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FIG. 2. The same as in Figure 1 for a germanium detector with threshold energy E t ^=12 keV. 
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FIG. 3. The same as in Figure 1 for a germanium detector with threshold energy Ej/j=4 keV 
and flat background b=0.1 cpd/kg/keV. 
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FIG. 4. The same as in Figure 1 for a TeO"2 detector with threshold energy E^=5 keV and flat 
background b=0.01 cpd/kg/keV. 
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FIG. 5. The same as in Figure 1 for an Nal detector with threshold energy Ej/j=2 keV and flat 
background b=0.1 cpd/kg/keV. 
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FIG. 6. Minimal exposure MTa as 



100 

M w (GeV) 

function of the WIMP 



1000 



mass mw required for the 



< 5 >=5.6 sensitivity plot to reach the exclusion plot of Ref. |5j to the WIMP-nucleon scalar 
cross section a^ n \ and calculated for a germanium detector with threshold energy E^=2 keV. 
Dotted, solid and dashed curves refer to v\ oc = 170, 220 and 270 km sec -1 respectively. For each 
value of v\ oc the curves calculated for a flat background b=0,0.01 and 0.1 cpd/kg/keV are given 
from bottom to top. 
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FIG. 8. The same as in Figure 6 for a Te02 detector with threshold energy E t ^=5 keV. 
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FIG. 9. The same as in Figure 6 for an Nal detector with threshold energy E t /j=2 keV. 
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